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ABSTRACT 



Context. Our knowledge of the optical spectra of Isolated Neutron Stars (INSs) is limited by their intrinsic faintness. 
Among the fourteen optically identified INSs, medium resolution spectra have been obtained only for a handful of 
objects. No spectrum has been published yet for the Vela pulsar (PSR B0833-45), the third brightest (V — 23.6) INS 
with an optical counterpart. Optical multi-band photometry underlines a flat continuum. 

Aims. In this work we present the first optical spectroscopy observations of the Vela pulsar, performed in the 4000- 
llOOOA spectral range. 

Methods. Our observations have been performed at the ESO Very Large Telescope (VLT) using the FORS2 instrument. 
Results. The spectrum of the Vela pulsar is characterized by a flat power- law F u oc v~ a with a — —0.04 zb 0.04 (4000- 
8000 A), which compares well with the values obtained from broad-band photometry This confirms, once more, that 
the optical emission of Vela is entirely of magnetospheric origin. 

Conclusions. The comparison between the optical spectral indeces of rot at ion- powered INSs does not show evidence 
for a spectral evolution suggesting that, as in the X-rays, the INS aging does not affect the spectral properties of the 
magnetospheric emission. At the same time, the optical spectral indeces are found to be nearly always flatter then the 
X-rays ones, clearly suggesting a general spectral turnover at lower energies. 
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1. Introduction 

The study of the optical emission properties of Isolated 
Neutron Stars (INSs) is hampered by their intrinsic faint- 
ness. For most of them the knowledge of the optical spec- 
trum is still based on the comparison between multi-band 
photometry measurements. We note here that this com- 
parison is affected by several uncertainties since photom- 
etry measurements are often compiled from the literature 
and, thus, are taken with different instruments and filters, 
calibrated with different photometric systems (167 used 
in total; see, e.g. Mono & Munari 200(0), indipendently 
corrected for the atmospheric and interstellar extinction 
and, last but not least, converted to spectral fluxes us- 
ing slightly different techniques. Only for five of the four- 
teen INSs with identified optical counterparts (see Mignani 
et al. 2004; Mignani 2005 for updated reviews) medium- 
resolution spectra have been obtained so far. The first one 
was the Crab pulsar (Oke 1969), the youngest (~ 1000 
years) and the brightest INS (V = 16.6), for which repeated 
spectroscopy observations have been performed (Nasuti et 
al. 1996; Sollerman et al. 2000). The second one was PSR 
B0540-69 (Hill et al. 1997, Serafimovich et al. 2004), the 
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youngest (~ 2000 years) and brightest (V = 22 A) INS after 
the Crab pulsar. Only after the advent of the generation of 
the 10 m class telescopes, it was possible to obtain spectra 
of the fainter (V ~ 25 - 25.6) and older (> 10 5 years) INSs 
such as Geminga (Martin et al. 1998) and PSR B0656+14 
(Zharikov et al. 2007), as well as of the radio-quiet INS RX 
J1856-3754 (van Kerkwijk & Kulkarni 2001). 

The best next target for optical spectroscopy is the 
Vela Pulsar (PSR B0833-45), the third optically brightest 
(V = 23.6) INS. Among the first radio pulsars discovered in 
the sixties, its optical identification was proposed by Lasker 
(1976) and confirmed by the detection of optical pulsations 
at the radio period (Wallace et al. 1977). The high pulsed 
fraction of the optical lightcurve (see also Gouiffes 1998), 
naturally pointed towards a non thermal, magnetospheric 
origin of the optical radiation. This was confirmed by multi- 
band photometry performed by Nasuti et al. (1997) with 
the NTT and by Mignani & Caraveo (2001) with the HST. 
The fluxes between - 3000 A and - 8000 A follow a flat- 
tish, power-law distribution (F v ex v~ a ;a = —0.2 ± 0.2). 
The non-thermal nature of the optical radiation of the Vela 
pulsar has been also confirmed by polarisation observa- 
tions performed with the VLT (Wagner & Seifert 2000; 
Mignani et al. 2007), pointing to a relatively high polar- 
isation level. The Vela pulsar has been recently observed 
in the IR with the VLT (Shibanov et al. 2003), and in the 
near-UV (Romani et al. 2005; Kargaltsev & Pavlov 2007) 
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with the HST. While its multi-band photometry now ex- 
tends from the IR to the near-UV, no optical spectrum 
has been published yet. Here, we report on the first spec- 
troscopy observations of the Vela pulsar performed with the 
VLT. Observations, data analysis and results are described 
in §2, while the interpretations are discussed in §3. 



2. Observations and Data Analysis 

2.1. Observations 

Spectroscopy observations of the Vela pulsar were per- 
formed in Service Mode in four different nights between 
December 2000 and January 2001 with the VLT at the 
ESO's Paranal Observatory (Chile). We have used FORS2 
(FOcal Reducer Spectrograph) , a multi-mode camera for 
imaging and long-slit/multi-object spectroscopy which is 
identical to FORS1 (Appenzeller et al. 1998) but with 
the CCD optmized to achieve a higher sensitivity in the 
Red part of the spectrum. At the time of the observations 
FORS2 was still mounting the original detector, a Tektronix 
2048x2046 pixels CCD with a plate scale of Of! 20 with 
the Standard Resolution (SR) collimator. The instrument 
was operated in its single-port, high-gain, read out mode 
which is the default one for long slit spectroscopy (LSS). 
To cover the wavelength interval 4000- 11 000 A the observa- 
tions were performed with two different grisms: the 300V 
(A A = 4500 - 8600 A) and the 300/ one (A A = 6000 - 1 1 100 
A), equipped with the order separation filter OG590. Both 
grisms have a dispersion of ~2.7A/pixel and a resolving 
power A/AA = 440. A slit width of 275 was used to collect 
as much flux as possible from the pulsar (see Fig. 1). 

Single LSS science exposures of 2800 s each were ob- 
tained for both grisms and repeated for cosmic rays fil- 
tering. A total of eight 2800 s exposures were taken for a 
total integration time of 11200 s for each grism. The com- 
plete log of the science observations is reported in Table 1. 
The first exposure sequence in the 300/ grism was aborted 
and is not considered in the following analysis. The see- 
ing conditions were always sub-arcsec (078 average) and 
the airmass below 1.2. The slit was oriented NS on the 
first night and always EW on the following nights. For 
each night, daytime calibrations (biases, darks, flatfields, 
arc lamp spectra) were taken to correct for instrumental 
effects and to perform the wavelength calibration. Multi 
object spectroscopy (MOS) observations of the spectropho- 
tometric standard stars LTT 3218, Feige 56 and Feige 67 
(Hamuy et al. 1994) were acquired at the beginning of each 
night (with the only exception of the third one) for flux cal- 
ibration. Because of their brightness, standard stars were 
observed in MOS mode (slitlet size 2270) to avoid flux losses 
which might occur when using the narrower LSS slit (2 / /5). 
Thanks to the sub-arcsec seeing conditions during the ob- 
servations it was possible to resolve the pulsar spectrum in 
each individual science exposure. As an example, we show 
in Fig. 2 one of the two 300V spectra taken on the first 
night, where the pulsar is clearly detected. Unfortunately, 
the pulsar is detected with a much lower significance in the 
300/ spectra due to the higher sky background and to the 
drop of the CCD sensitivity towards longer wavelengths 
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Fig. 1. VLT FORS2 raw £?-band acquisition image of the 
Vela pulsar. North to the top, East to the left. The 2"5 LSS 
slit is overplotted, with the orientation NS as used in the 
first night (see Table 1). The position of the Vela pulsar 
(barely detected in the 60 s exposure acquisition image) 
is marked by the arrow. The difference between the four 
quadrants corresponds to the different CCD readout ports. 

Table 1. Summary of the spectroscopy observations of the 
Vela pulsar taken with FORS2. The columns give the ob- 
serving dates, the used grism, the exposure times, the num- 
ber of exposures per pointing N, and the average seeing and 
airmass values during each exposure sequence. 



yyyy-mm-dd 


Grism 


Time (s) 


N 


Seeing 


Airmass 


2000-12-03 


300V 


2800 


2 


or 79 


1.14 


2001-01-23 


300V 


2800 


2 


0760 


1.07 


2001-02-25 


3001 


1400 


1 


0789 


1.09 


2001-01-30 


3001 


2800 


4 


0773 


1.11 



2.2. Data Analysis 

Spectral data reduction and calibration were performed us- 
ing standard suites of tools for treating CCD data available 
in the MIDAS and IRAF sofware packages. Master cali- 
bration frames (master biases and darks, flux-normalized 
fiat fields) were provided by the FORS data reduction 
pipelined Science and standard stars spectra, as well as 
arc lamp spectra, were then bias and dark-subtracted, and 
flatfielded. Wavelength calibration was computed from the 
reduced arc lamp spectra by fitting a second order polyno- 
mial, yielding an rms of 0.4 A/pixel. The wavelength cali- 
bration was then applied row-by-row to both the standard 
stars and science spectra. For each grism, the flux normal- 
ization was computed from the extracted one-dimensional 
standard star spectrum using the available flux reference ta- 
ble and applying the atmospheric extinction correction with 
the extinction curves measured for the Paranal Observatory 
(Patat 2004). For each grism, single science spectra were 
finally co-added and cosmic rays removed. The pulsar one- 
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dimensional spectrum was then extracted from each of the 
co-added science spectra using a 4 pixel wide (Of! 8) ex- 
traction window centered on the pulsar position where the 
S/N is higher. The sky background was computed and sub- 
tracted from two 6 pixel wide (1//2) regions immediately 
adjacent to the pulsar spectrum. For each extracted one- 
dimensional spectrum of the pulsar we finally applied a flux 
correction factor, computed from the extracted spectra of 
brighter field stars, to account for the width difference be- 
tween the extraction window and the object's PSF. The 
spectra were corrected for the interstellar extinction using 
as a reference E(B — V) = 0.05 obtained from the spectral 
fit to the XMM X-ray spectrum (e.g. Manzali et al. 2007). 
The two one-dimensional spectra of the pulsar were finally 
merged to yield a single spectrum across the 4000- 11000 A 
wavelength range. 

2.3. Results 

The final Vela pulsar's spectrum is shown in Fig. 3. Due 
to the different detection significance of the pulsar in the 
two grisms, the S/N of the resulting (unbinned) spectrum is 
rather unhomogeneous. In particular, we estimated a S/N 
~ 5 (per pixel) in the 4000-5500A wavelength range and 
a S/N as low as ~ 1 near/above 8000A. At longer wave- 
lengths the spectrum is entirely dominated by noise. For 
this reason in the following analyses we have considered 
only the part of the spectrum between 4000 A and 8000 
A. Clearly, the spectrum consists of a pure flat contin- 
uum, with no evidence for emission or absorption feau- 
tures. For comparison, we have overplotted the spectral 
fluxes recomputed from the optical multi-band photome- 
try (Nasuti et al. 1997; Mignani & Caraveo 2001), obtained 
between ~ 4000 A and ~ 8000 A. To be consistent, we 
have re-corrected these fluxes for the same value of the in- 
terstellar extinction used in this work (Ay =0.16), while 
both Nasuti et al. (1997) and Mignani & Caraveo (2001) 
have used Ay=0A. For the interstellar extinction correc- 
tion we have used the coefficients of Fitzpatrick (1999). As 
seen, within the cross-calibration uncertainties, the multi- 
band photometry and the spectroscopy flux measurements 
are substantially consistent with each other. A putative dip 
at 6500 A was hinted in the multi-band photometry data 
of Mignani & Caraveo (2001). However, its existence was 
not confirmed by the VLT photometry of Shibanov et al. 
(2003). No evidence for such a dip is found in our spectral 
data either. We thus conclude that it was just the result of 
a cross-calibration problem in the multi-band photometry 
of Mignani & Caraveo (2001). 

We have fitted the FORS2 spectrum between 4000 and 
8000 A with a power-law (F v oc u~ a ) and we derived 



a spectral index a = —0.04 ± 0.04. This value super- 
seeds the one obtained from the multi-band photometry 
of Mignani & Caraveo (2001) on the base of only five spec- 
tral flux measurements. We note that our spectral index 
is somewhat flatter than the one obtained by Shibanov 
et al. (2003) by fitting the IR-to-optical spectral fluxes 
(a = 0.12 ± 0.05), while it compares better with the spec- 
tral index a = 0.01=b0.02, obtained by Kargaltsev & Pavlov 
(2007) by extending the fit to the near-UV. Interestingly, 
the optical power-law is below the extrapolation of the 
power-law component (a = 1.2 ± 0.3) used to fit the XMM 
X-ray spectrum (Manzali et al. 2007; see also Kargaltsev 
& Pavlov 2007), which clearly suggests that the magne- 
tospheric emission of the Vela pulsar features a spectral 
turnover at longer wavelengths. 



3. Discussion 

We have compared the spectrum of the Vela pulsar with 
those of the other rotation-powered INSs for which either 
medium-resolution spectroscopy or multi-band photome- 
try is available (Fig. 4). Apparently, the complexity of the 
spectral flux distribution grows with the INS age. For the 
young objects the optical spectral energy distribution is 
dominated by a flat power-law continuum which brings the 
signature of synchrotron radiation produced by relativis- 
ts charged particles in the neutron star's magnetosphere. 
For the older ones, a Ray leigh- Jeans component, ascribed 
to thermal radiation from the cooling neutron star's sur- 
face, is also present. In all cases, no statistically signifi- 
cant evidence of emission or absorption features is found. 
For all the objects in Fig. 4, Tab. 2 reports the measured 
spectral index ao of the power-law component, either de- 
rived from spectroscopy or from broad-band photometry. In 
those cases where both values are available the agreement 
is rather good, with the only exception of PSR B0540-69. 
However, as discussed in Serafimovich et al. (2004), its spec- 
trum is unrecoverably polluted by the background of the 
surrounding, bright compact (~ 4") synchrotron nebula. 
The comparison of the optical spectral indeces shows that 
it is difficult to find clear spectral templates for different 
INS groupfl This is shown in Fig. 5 (top left panel), where 
we have plotted the optical spectral index ao as a func- 
tion of the INS spin-down age. When available, we have 
assumed as a reference the spectral index obtained from 
spectroscopy, with the obvious exception of PSR B0540-69 
(see above). 



4 We warn here that the values of ao have been computed 
over slightly different wavelength ranges, which makes a direct 
comparison more uncertain 
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Fig. 3. Extracted one-dimensional spectrum of the Vela pulsar obtained by merging the 300V and 300/ spectra (shown 
in green and cyan, respectively) after wavelength/flux calibration, extinction correction, and coaddition. A rebinning of 
a factor 3 has been applied to both spectra. The spikes are due to defects in the subtraction of bright sky lines and of 
residual cosmic ray hits. The red points correspond to broad-band photometry (Nasuti et al. 1997; Mignani & Caraveo 
2001). The red dotted line represents the best fitting power-law (a = —0.04 ± 0.04) to the continuum. Since at longer 
wavelengths the spectrum is dominated by the sky background, the fit has been performed only between 4000 and 8000 
A. 



Table 2. Summary of the measured optical/IR spectral index ao of rotation-powered INSs as obtained from spectroscopy 
(column three) and multi-band photometry (column five). Columns 2-4 give (in logharitmic units) the pulsar spin down 
age r, the magnetic field B, and the rotational energy loss E. The last column gives the phase- averaged X-ray spectral 
index ax- 



Name 


Log(r) 
yrs 


Log(B) 
(G) 


Log(£) 
(erg s x ) 




AA 

(A) 


&o, P h 


AA 

(A) 


ax 


Ref. 


Crab 


3.1 


12.58 


38.65 


-0.11 ±0.04 


3300-9250 






1.079 ±0.004 


1,2 


B1509-58 


3.2 


13.19 


37.25 






0.5 


6000-18000 


0.4 ±0.5 


3,4 


B0540-69 


3.2 


12.70 


38.17 


1.88 ±0.01 


3605-6060 


1.07 ±0.2 


3300-8000 


0.88 ±0.11 


5,4 


Vela 


4.1 


12.53 


36.84 


-0.04 ±0.04 


4000-8000 


0.01 ±0.02 


1500-16000 


1.2 ±0.3 


6,7,8 


B0656+14 


5.0 


12.67 


34.58 


0.2 ±0.2 


4600-7000 


0.41 


4300-18000 


1.1 ±0.3 


9,10,11 


Geminga 


5.5 


12.21 


34.51 


0.8 ±0.5 


3700-8000 


0.46 


4300-16000 


0.7 ±0.1 


12,10,11 


B1929+10 


6.5 


11.71 


33.59 






0.5 ±0.5 


1700-3400 


1.7 ±0.1 


13,14 


B0950+08 


7.2 


11.39 


32.75 






0.65 ±0.4 


3600-8000 


0.9 ±0.1 


15,16 



(I) Sollerman et al. (2000), (2) Kirsch et al. (2006), (3) Kaplan & Moon (2006), (4) Gotthelf (2003), (5) Serafimovich et al. (2004), 
(6) this work, (7) Kargaltsev & Pavlov (2007), (8) Manzali et al. (2007), (9) Zharikov et al.(2007), (10) Koptsevich et al. (2001), 

(II) De Luca et al. (2005), (12) Martin et al. (1998), (13) Mignani et al. (2002), (14) Becker et al. (2006), (15) Zharikov et al. 
(2004), (16) Becker et al. (2004) 



Although the spectral index value seems generally to 
correlate with the spin-down age, with ao ~ (0.12 ± 
0.04)Log(r), the errors are such that most of the points 
are also consistent with a costant distribution around the 
average value < ao >= 0.44 ± 0.4. Thus, there is no firm 
evidence for an evolution of the pulsars' optical spectral in- 
dex over four age decades. This is similar to what is found 
in the X-rays, where the available data also suggest that the 
value of the spectral index ax does not depend on the spin- 
down age (see, e.g. Becker & Triimper 1997). A possible anti 
correlation between the optical spectral index and the ro- 
tational energy loss E, with ao ~ (-0.07 ± 0.02)Log(E), 
can also be recognized in the data (Fig. 5, top right panel). 
However, as before, the large errors do not allow to draw 
any firm conclusion. Finally, no trend can be recognized 



between the optical spectral index and the magnetic field 
B (Fig. 5, bottom left panel). 

As it has been pointed out (e.g. Mignani et al. 2004; 
Serafimovich et al. 2004), the optical and X-ray magne- 
tospheric emission of rotation-powered pulsars are almost 
never described by the same spectral parameters. This is 
clearly shown in Fig. 5 (bottom right panel), where we have 
plotted the optical spectral index ao vs. the X-ray one ax- 
As seen, the optical power-law index is often less steep than 
the X-ray one, suggesting that the spectra of rotation-power 
pulsars undergo a general turnover at lower energies. In par- 
ticular, for the Crab the slope of the power-law swaps from 
positive to negative when passing from the X-ray to the op- 
tical/IR domain, underlying an even more marked spectral 
turnover. Furthermore, no general positive or negative cor- 
relation is found between ao and olx, which indicates that 
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Fig. 5. Values of the optical spectral index ao plotted as a function of the INS spin-down age, rotational energy loss 
E, the magnetic field B, and the X-ray spectral index ax (top left to bottom right). The red dashed lines (top panels) 
represent the linear fit to the points. The dashed line (bottom right panel) marks the case ao = olx- In all panels the 
dot-dashed and dotted lines correspond to the average value (< ao >= 0.44 ± 0.4) of the optical spectral index and to 
its la variation, respectively. 



the X-ray and optical magnetospheric emission, although 
likely produced by the same physical processes, are not di- 
rectly correlated. This finding is particularly interesting in 
comparison with the strong correlation between the optical 
and X-ray luminosities of rotation-powered pulsars noticed 
by Zharikov et al. (2004) and by Zavlin and Pavlov (2004). 

4. Summary 

We have presented the first optical spectroscopy observa- 
tions of the Vela pulsar. The pulsar's spectrum is char- 
acterized by a flat power-law with spectral index a = 
—0.04 ±0.04, consistent with the values derived from broad- 
band photometry (Mignani & Caraveo 2001; Shibanov et al. 
2003; Kargaltsev & Pavlov 2007). We have compared the 
newly derived optical spectral index of Vela with those of 
all rotation-powered INSs for which a power-law compo- 



nent has been identified in the optical/IR spectrum. While 
a trend can be recognized in the data, the large errors on 
the spectral index values for most of the objects prevent any 
claim for an evolution of the magnetospheric emission prop- 
erties over four age decades. We also found no evidence for a 
correlation between the optical and X-ray spectral indeces, 
which indicates that the X-ray and optical magnetospheric 
emissions are not directly correlated. However, we showed 
that in the majority of cases, the optical spectral index is 
flatter than the X-ray one, suggesting a spectral turnover 
in the INSs' spectra at low energies. 
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Fig. 4. Spectral flux distribution of all rotation-powered 
pulsars for which either medium-resolution spectroscopy 
or multi-band photometry is available (see Table 2). From 
top to bottom, objects are sorted according to increasing 
spin-down age. This figure has been updated from Fig. 1 of 
Zharikov et al. (2007), by including the IR fluxes of PSR 
B1509-58 (Kaplan & Moon 2006), the Vela spectrum (this 
work), and the near-UV fluxes of Vela (Romani et al. 2005), 
Geminga (Kargaltsev et al. 2005) and PSR B0656+14 
(Shibanov et al. 2005; Kargaltsev & Pavlov 2007). 
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